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Abstract. We argue that the collinear factorization of the fragmentation 
functions in high energy nuclear collisions breaks down at transverse momenta 
Pt ^ Qs/g due to high parton densities in the colliding hadrons and/or nuclei. 
We find that gluon recombination dominates in that pt region. We calculate the 
inclusive cross-section for it meson and nucleon production using the low energy 
theorems for the scale anomaly in QCD, and compare our quantitative baryon- 
to-meson ratio to the RHIC data. 



PACS numbers: 13.85.Ni, 24.85. +p, 25.75.-q 

1. Introduction 

It is believed that the strong interactions at high energies/densities are described by 
the Colour Glass Condensate (CGC), which provides the theoretical framework for 
the description of the saturation of the gluon fields in the hadronic and nuclear wave 
functions pQ. The transition to the saturation region is marked by the saturation 
scale Q s , which is the two-dimensional parton colour charge density in the infinite 
momentum frame of the hadron wave function. The typical momentum scale of the 
partons in the saturation regime is in the order of Q s . 

We argued in [2] that the process of gluon fusion in high parton density medium 
should dominate in the J PC = ++ singlet state production channel at pt <C Q s / ' ^/al 
region. Contribution of this channel to the total multiplicity is proportional to the 
spectral density pg(M) of the correlator (0\T{9^(x), 6^(0)}|0). In the perturbation 
theory, valid at large invariant masses, the lowest order contribution to the correlator 
arises from the two gluon amplitude | (gg\9^(x)\0) | 2 . On the other hand, at low 
invariant masses, the main contribution comes from the excitation of hadrons through 
the \(hh- ■ • \0fi(x)\0}\ 2 amplitude. We can then apply anomaly matching to infer the 
hadron production from the two-gluon production amplitude in the same channel by 
weighting it with an appropriate weight factor, which can be calculated in the chiral 
perturbation theory. 

2. Gluon pair production in the scale anomaly channel 

The hadron production model discussed in this paper consists of two stages: (i) 
the production of a gluon pair in J PC = ++ and colour singlet state, and 
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Figure 1. Production of gluon and pion pairs in the J PC = ++ channel. 



(ii) recombination of this pair into final state hadrons via the anomaly matching 
mechanism [3J. It is shown in [2] that this gluon fusion process exceeds the collinear 
factorization channel and dominates in the double inclusive gluon production. The 
relevant diagrams are given in figure [1] In the approximation that gluons are typically 
produced back-to-back, the spectrum of a single gluon with transverse momentum px 
is given by [4] 

daf + n 4 f ip T , , ^ , 2N [ 2 ^dM 2 p qq (M) 

where <Pa(x, q 2 ) is the unintegrated gluon distribution function and p gg is the spectral 
density in the lowest order perturbation theory [5] 
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3. Hadronic contributions to the spectral density 

At small invariant masses of the produced gluons, the leading contribution arises from 
the hadronic excitations and can be calculated in the chiral perturbation theory. The 
chiral Lagrangian for pions is given by 

C n = LTrdpUd»Ut + ±mlfiTr(U + tf), (3) 

where U = exp{2i7r//7r}, it = ir a T a and T a are the SU(2) generators normalized by 
Tr (T a T b ) — \5 ah ■ The non-perturbative contribution for pions to the spectral density 
can be calculated as p 7r ^{M 2 ) = ^-^M 4 , where M is the invariant mass of the initial 
gluon pair [5J. For nucleons, the chiral Lagrangian reads 

CN = m N tptp, (4) 

which gives the spectral density for nucleons pjqfj (M 2 ) = ^ m 2 N M 2 . The spectral 
density for nucleon kinetic term has a much weaker dependence on the invariant mass 
M, and cannot compete with the M 4 dependence for pion kinetic term. Therefore, 
we need to include the pion-nucleon interaction contribution to the spectral density, 
which manifest itself as the following Yukawa Lagrangian 5 

2im N g A ^- 

Li = 77^75 1 N . (5) 

Jrr 



Baryon enhancement in high- density QCD and relativistic heavy ion collisions 3 



The spectral density for pion-nucleon interactions scales as M 4 for high invariant 
masses by evaluating the three-particle final-state diagram [BJ, 
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which ensures that nucleons can be produced with the same orders of magnitude 
as pions. Their ratio of total cross sections will approach an asymptotic value for 
high enough colliding energy. Hadron production cross sections can be obtained by 
replacing the spectral density for double-gluon production in ([T]) by the corresponding 
ones for hadrons in the chiral perturbation theory. 



4. Baryon-to-meson ratio in AuAu collisions at RHIC 

We will now apply our formulas obtained in the previous sections to calculate the p 
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Here the single hadron production cross section d °^± is added since it is expected 
that hard processes will dominate particle productions at high transverse momentum. 
The coefficients of the soft production channel suggest an isospin symmetry 

The single gluon production cross section decreases as l/p? and will dominate, 
at some high pt, the one in JT]), which drops as \jp%- The hadron cross section at 
high pt is then given by the convolution of gluon production cross section with an 
appropriate fragmentation function [T\ 

do { h ] _ f dz da g (p T /z) g/h 



dp 2 T dy~J z 2 dp\dy F ^ Z ' Pt) - (8) 

There is now only one parameter to be specified which is the saturation scale Q s . 
According to the KLN model [21 [8] it is given to be the sum of the saturation scales 
of each nuclei, which at RHIC is Q 2 ps 2 x 0.284 N^ t GeV 2 at midrapidity. 

The baryon-to-meson ratio (BMR) is plotted as a function of the hadron 
transverse momentum and the number of participants at ^/snn = 200 GeV in figure 
We can see the enhancement of the ratio in both cases. Since baryon production has 
a higher threshold energy, the available phase space is much smaller than for mesons. 
As a result it is suppressed at small pr and 7V par t compared to meson production. 
By going from peripheral to central collisions, the saturation scale increases. Hence 
the increase of a typical gluon transverse momentum makes possible the production 
of baryons. At px ^> Q s our recombination production mechanism gives way to the 
traditional gluon fragmentation in accordance with the parton-hadron duality, and 
the baryon-to-meson ratio is now mainly determined by the pQCD calculations. 

Our calculations showed that the cold nuclear effects due to the gluon saturation 
are reasonably described in AuAu collisions. We believe that this mechanism also 
contributes to the fast chemical equilibration in the produced dense medium in AuAu 
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Figure 2. Baryon-to-meson ratio in AuAu collisions at ^/snn = 200 GeV as 
a function of transverse momentum (a) and centrality (b). (a) Data are 0-12% 
central AuAu collisions [9]. Two lines correspond to the following relationships 
between the relative momentum k and the hadron momentum px- Dashed line: 
k 2 4p|. Solid line: k 2 RJ 4(p T - A) 2 , with A = 0.25 GeV being a small 
offset momentum, (b) Data are PHENIX preliminary results for integrated 3 
GeV < p T < 4 GeV region. 



collisions by significantly enhancing the baryon yield as compared to the pQCD 
calculations. We suggest that the ratio is close to the thermal equilibrium value 
already at the initial time. The final state interaction in the hot medium will achieve 
a thermal equilibration and may further increase the baryon-to-meson ratio. 

5. Summary 

We have calculated the baryon-to-meson ratio due to cold nuclear effects in the high 
parton density QCD. We found that a gluon fusion mechanism dominates hadron 
production and it will significantly enhance the BMR already in the initial time. 
Combined with pQCD calculation for high transverse momentum region, our model 
can reasonably describe part of the RHIC AuAu data. Of course, the thermal 
equilibration is achieved only through the final state interactions. 

Acknowledgments 

I would like to thank my collaborators Kirill Tuchin and Dmitri Kharzeev. 
References 

[1] Iancu E and Venugopalan R 2003 Quark Gluon Plasma 3 ed Hwa R C and Wang X N (Singapore: 
World Scientific) p249 
Jalilian-Marian J and Kovchegov Y V 2006 Prog. Part. Nucl. Phys. 56 104 
[2] Li Y and Tuchin K Phys. Rev. D 75 074022 
[3] Kharzeev D and Levin E Nucl. Phys. B 578 351 
[4] Kharzeev D, Li Y and Tuchin K, to be submitted 

[5] Weinberg S 1996 The Quantum Theory of Fields: Modern Applications vol 2 (Cambridge: 

Cambridge University Press) chapter 19 
[6] Yao W-M 2006 J. Phys. G: Nucl. Part. Phys. 33 1 
[7] Kharzeev D, Kovchegov Y V and Tuchin K 2004 Phys. Lett. B 599 23 
[8] Kharzeev D and Nardi M 2001 Phys. Lett. B 507 121 
[9] Abelev B I est al (STAR Collaboration) 2006 Phys. Rev. Lett. 97 152301 



